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SZ effect from Dark Matter annihilation 
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Abstract. We derive in this Letter the SZ effect induced by the secondary electrons produced in the annihilation 
of Weakly Interacting Massive Particles (assumed here to be neutralinos) in gravitationally bound structures 
dominated by Cold Dark Matter (CDM). We show that the DM induced SZ effect has a specific spectral shape 
and an amplitude which increases for decreasing neutralino mass M x . The available SZ data on the Coma cluster 
set an upper limit on the quantity {aV)An x which can be combined with the WMAP constraints on Q m h 2 to 
restrict the available neutralino models in the (oV)a — M x plane. We delineate various potential applications of 
this method to constraint the physical properties of the Dark Matter particles from the study of galaxy clusters 
and dwarf spheroidal galaxies in the light of the next coming high-sensitivity SZ experiments. 
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1. Introduction 

Dark Matter (DM) annihilation in the halo of galaxies 
and galaxy clusters have relevant astrophysical implica- 
tions. In fact, if DM is constituted by weakly interacting 
massive particles (WIMPs), their annihilation can pro- 
duce direct and indirect signals such as observable fluxes 
of positrons (e.g., Silk & Srednicki 1984, Kamionkowski 
& Turner 1991, Baltz & Edsjo 1999) antiprotons (e.g., 
Bottino ct al. 1998), gamma rays (e.g., Chardonnet et 
al. 1995, Colafrancesco & Mele 2001), neutrinos (e.g., 
Gondolo & Silk 1999, Hooper & Silk 2004), radio emis- 
sion (e.g., Colafrancesco & Mele 2001), heating of the hot 
intra-cluster gas (Totani 2004, Colafrancesco 2004). We do 
not have yet however, at present, a definite detection of 
these emission features originating from DM annihilation. 

As an alternative strategy, we explore here the conse- 
quences of the Compton scattering between the secondary 
electrons produced from the WIMP anihilation in massive 
DM halos, like galaxy clusters, and the CMB photon field. 

Galaxy clusters are gravitationally dominated by Cold 
Dark Matter for which the leading candidate is the lightest 
supersymmetric (SUSY) particle, plausibly the neutralino 
X- Experimental and theoretical considerations for having 
a cosmologically relevant neutralino DM lead to bound its 
mass M x in the range between a few GeV to a few hun- 
dreds of GeV (e.g., Bottino et al. 2003, Belanger 2003). 
The decays of neutralino annihilation products (fermions, 
bosons, etc.) yield, among other particles, energetic elec- 
trons and positrons up to energies comparable to the neu- 
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tralino mass. Here we notice that these energetic electrons 
and positrons (hereafter we will refer to these particles 
as electrons because their distinction is not essential for 
our purpouses) can interact with the CMB photons and 
up-scatter them to higher frequencies producing a pecu- 
liar Sunyaev-ZePdovich (1972, 1980, hereafter SZ) effect 
with specific spectral and spatial features. In this Letter, 
we will describe the specific feature of the SZ effect pro- 
duced by DM annihilation, SZ£>a/i m galaxy clusters and 
we will discuss the possibility to disentangle such specific 
SZjjjvf effect from the other sources of SZ effect which 
are present in the same structures. We will finally discuss 
the future experimental prospects for the detection of the 
SZdm effect in DM halos. The relevant physical quanti- 
ties are calculated using Hq = 70 km s _1 Mpc^ 1 and a 
flat, vacuum-dominated CDM (fi m = 0.3, ft\ = 0.7) cos- 
mological model. 

2. The SZ effect from DM annihilation 

The calculation of the secondary electron spectrum from 
XX annihilation in galaxy clusters has been already pre- 
sented in details by Colafrancesco & Mele (2001) and here 
we will only recall the relevant steps which are necessary 
for the present purpouses. 

2.1. Neutralino annihilation in galaxy clusters 

Following Colafrancesco & Mele (2001) we assume a spher- 
ical, uniform top-hat model for the collapse of a cluster 
with an average DM density p — A(f2o, z)pb(l — f g ) where 
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the background density is pi, = Ct m p c and A(£lo, z) is the 
non-linear density contrast of the halo which virializes at 
rcdshift z. Here f g is the gas mass fraction. Redistributing 
the total mass, M = 4^ pi? 3 , found within the radius 
R p = pr c (expressed in terms of the cluster core radius 
r c ), according to a density profile p(r) = pag{r), we find 

- R 3 

that the central total density is po — ffpfj , where 

I(Rp) — Jq P drr 2 g(r). In our phenomenological approach, 
we consider mainly the case a central cusp profile, as in- 
dicated by the available N-body simulations. The central 
cusp model described by a profile: 

ri-t 



9(r) = 



1 + 



(1) 



where I(p, 77, £) = 
3, the central 



gives a central density p = § j^^j 
J£dxx 2 - r '(l + x)^. With r, = 1 and £ 
cusp model corresponds to the "universal density profile" 
which Navarro, Frenk & White (1997) showed to be a good 
description of cluster DM halos in N-body simulations of 
CDM hierarchical clustering. Assuming that the DM den- 
sity scales like the total cluster density, the central DM 
number density is: 

„3 



n x,o = 



np 
3T 



where 



■ 4.2 • 10 _b cm 



- 3 n x h 2 



M x 


-1 


rA(fio,*)i 


100 GeV_ 




400 



(2) 



(3) 



The x annihilation rate in a DM halo is R — n x (r)(aV) a , 
where n x (r) — n x ^g(r) is the neutralino number density 
and ((jV)a is the XX annihilation cross section averaged 
over a thermal velocity distribution at freeze-out tempera- 
ture. Although the XX annihilation cross section is a non- 
trivial function of the mass and physical composition of 
the neutralino, to our purpose it suffices to recall that the 
X relic density is approximately given by: 

-1 



n x h 2 



3 X 10 



~ 27 cm 3 s 



,VU <4) 

(Jungman ct al. 1996). Hence, for the XX annihilation, 
we shall assume a total cross section of {<jV)a ~ 2.6 • 
10~ 26 cm 3 s _1 to be consistent with the value Vl m h 2 ~ 
0.116 derived from WMAP (Bennett ct al. 2003). For val- 
ues of £l m h 2 in the range 0.085 — 0.152, the annihilation 
cross section is fixed to within a factor less than three. 
Detailed studies of the relic neutralino annihilation (Edsjo 
1997) show that the above value is well inside the allowed 
range predicted in supersymmetric theories for a wide 
choice of masses and physical compositions of neutralinos 
that can be relevant as CDM candidates. Enhancing (sup- 
pressing) the x annihilation rate will have on our results 
the simple effect of rescaling the final electron spectra by 
the same enhancement (suppression) factor. 

Neutralinos which annihilate inside a DM halo produce 
quarks, leptons, vector bosons and Higgs bosons, depend- 
ing on their mass and physical composition. Electrons are 



then produced from the decay of the final heavy fcrmions 
and bosons (monochromatic electrons, with energy about 
M x , coming from the direct channel XX ~ * ee i are i n gen- 
eral much suppressed, Turner & Wilczek 1990). The dif- 
ferent composition of the XX annihilation final state will 
in general affect the form of the final electron spectrum. 
Analytical expressions for the e ± spectrum has been given 
by Colafrancesco & Mele (2001) and we refer to this paper 
for further details. 

The time evolution of the electron spectrum is given 
by the transport equation: 



dn e (E, r) 
dt 



d_ 

dE 



n e (E,r)b(E) 



Qe(E,r) 



(5) 



where n e {E,r) is the equilibrium spectrum at distance 
r from the cluster center for the electrons with en- 
ergy E. The source electron spectrum, Q e (E,r), rapidly 
reaches its equilibrium configuration mainly due to syn- 
chrotron and Inverse Compton Scattering losses at ener- 
gies E <; 150 MeV and to Coulomb losses at smaller ener- 
gies (e.g., Colafrancesco & Mele 2001). Since these energy 
losses are efficient in the cluster atmosphere and DM an- 
nihilation continuously refills the electron spectrum, the 
population of high-energy electrons can be described by a 
stationary transport equation 



d_ 

dE 



n e (E,r)b e (E) 



Qe(E,r) 



(6) 



from which the equilibrium spectrum can be calculated. 
Here, the function b e (E) gives the energy loss per unit 
time at energy E 



b e (E) = b (B, 



E 
GeV 



bco 



(7) 



where &o(£ M ) = (2.5 • 10~ 17 + 2.54 • 10~ 18 .B 2 ) and b Co ui = 
7 x 10~ 16 [n(r)/lcm~ 3 ], if b e is given in units of GeV/s. 
The source spectra we derived (see Fig.l in Colafrancesco 
& Mele 2001) arise from an analytic approximation of the 
exact shape of the electron spectrum that tries to cope 
with the details of the quarks and leptons decays and of 
the hadronization of the decay products. Detailed elec- 
tron spectra can also be obtained by using state-of-the-art 
Monte Carlo simulations, although the analytical approxi- 
mations used here can resume the relevant aspects of more 
detailed studies. 

2.2. The SZ effect from x~ X annihilation products 

The generalized expression for the SZ effect which is valid 
in the Thomson limit for a generic electron population in 
the relativistic limit and includes also the effects of mul- 
tiple scatterings and the combination with other electron 
population in the cluster atmospheres has been derived by 
Colafrancesco et al. (2003). This approach is the one that 
should be properly used to calculate the specific SZdm 
effect induced by the secondary electrons produced by XX 
annihilation. Here we do not repeat the description of the 
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analytical technique and we refer to the general analy- 
sis described in Colafrancesco et al. (2003). According to 
these results, the DM induced spectral distortion can be 
written as 



at i \ A^bTo) . 
A/dm N = 2 y DM g(x) 

(he)'' 



Vdm 



PdmM , 



<?0) 



m e c 
(k B T e ) 



1 



io(xe s )P(s)ds — io(x) 



in terms of the photon redistribution function P(s) and 
of i (x) = 2(k B T ) 3 /(he) 2 ■ x 3 /(e x - 1), where we defined 
the quantity 



i 

dpf e {p)^pv{p)m e c 



0.0006 



(8) 



where To is the CMB temperature and the 
Comptonization parameter yr>M is given by 



0.0002 



(9) 



0.0002 



in terms of the pressure Pdm contributed by the sec- 
ondary electrons produced by neutralino annihilation. The 
quantity udm oc (aV)An 2 , and scales as oc (<jV)aM~ 2 , 
providing an increasing pressure Pdm and optical depth 
Tdm — a T J d£n e for decreasing values of the neutralino 
mass M x . The function g(x), with x = hv/ksTfj, for the 
DM produced secondary electron population can be writ- 
ten as 




(10) 



(see Colafrancesco et al. 2003) which is the analogous 
of the average temperature for a thermal population 
(for a thermal electron distribution (ksT e ) = ksT e ob- 
tains, in fact). The photon redistribution function P(s) = 
j dpf e (p)P s (s;p) with s = ln{y' jv), in terms of the CMB 
photon frequency increase factor v' jv = ^j 2 — |, depends 
on the electron momentum (p) distribution, / e (p), pro- 
duced by XX annihilation. 

Figffl shows a comparison between the SZ effect ob- 
tained for a thermal electron distribution with k B T = 8.2 
keV and that evaluated for the XX annihilation for dif- 
ferent values of M x . A major difference between the two 
spectral functions is the different position of the zero of 
the SZ effect which is moved to higher frequencies in the 
case of the DM produced electrons with respect to the case 
of the thermal distribution. As a consequence, the SZdm 
effect appears as a negative contribution to the overall SZ 
effect at all the frequencies which are relevant for the SZ 
experiments, x ~ 0.5 — 10. Note that the amplitude of the 
SZ/)a/ effect increases for decreasing values of the neu- 
tralino mass since the quantity yoM °c (&V)An x ~ M~ 2 
which accounts for the scaling n x <~ M" 1 . A further M~ x 
dependence is present in the electron source function from 
XX annihilation (see, e.g., Colafrancesco & Mele 2001, 
their appendix). 



Fig. 1. The overall SZ effect in Coma produced by the com- 
bination of various electron populations: thermal hot gas with 
k B T = 8.2 keV and r = 4.9 • 1CT 3 (solid blue curve) which 
best fits the available SZ data (DePetris et al. 2003); relativis- 
tic electrons which best fit the radio-halo spectrum (yellow 
curve) provide a small additional SZ effect (Colafrancesco et 
al. 2003); warm gas with k B T « 0.1 keV and n « 10" 3 cm" 3 
(cyan curve) provides a small SZ effect due to its low pressure 
(Colafrancesco 2004); DM produced secondary electrons with 
M x = 10 (black dotted curve), 20 GeV (black solid curve) and 
30 GeV (dashed solid curve) . A pure-gaugino \ reference model 
is assumed in the computations. The relative overall SZ effect is 
(Hi own as the dotted, solid and dashed red curves, respectively. 
A zero peculiar velocity of Coma is assumed consistently with 
the available limits (Bernardi et al. 2002). SZ data are from 
OVRO (magenta), WMAP (cyan) and M1TO (blue). 



2.3. Application to galaxy clusters 

The SZ effect in galaxy clusters is dominated by the 
Compton scattering produced by the hot, thermal intra- 
cluster gas. However, it has been shown that there are 
additional sources of SZ effect which might be produced 
by other electron populations of thermal (warm IC gas, 
see e.g. Colafrancesco 2004) and non-thermal (relativistic 
electrons producing, e.g., the cluster radio-halo emission, 
Colafrancesco et al. 2003, Colafrancesco 2004) origin resid- 
ing in the cluster atmosphere. The SZ effect unavoidably 
induced by DM annihilation in the same cluster combines 
with the previous SZ effects to produce an overall SZ ef- 
fect which, hence, retains a complex information about 
the physical structure of the cluster atmosphere. We will 
refer, in the following, to the case of the Coma cluster for 
which there are several SZ observations available at differ- 
ent frequencies, from ~ 30 GHz to ~ 350 GHz. The Coma 
cluster atmosphere is a complex environment which hosts 
various electronic populations: a hot, thermal gas with 
fcfjT w 8.2 keV producing the bulk of the X-ray emis- 
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100 1000 

Log W x (GeV) 

Fig. 2. The constraints on the ((tV)a - M x plane set by the 
SZ effect from Coma. The heavily dashed area is excluded by 
the analysis of the SZdm in Coma. The SZ constraints are 
combined with the WMAP constraint 0.085 < Q. m h 2 < 0.152 
(Bennett et al. 2003) which are translated on the quantity 
{oV)a through eq.l0J. 



sion in the 1-10 keV region; a population of relativistic 
electrons with a non-thermal spectrum, n re i oc E~ x with 
x ~ 3.5 — 4, and energy E e 



16.4GeVB / ; 1/2 (i//GHz) 1 /2 ; 



which are responsible for the radio-halo synchrotron emis- 
sion (here is the magnetic field in units of /zG, see, e.g., 
Colafrancesco 2004 for a review) ; a possible warm gas halo 
with ksT w 0.1 keV claimed to be responsible for the ob- 
served soft X-ray/EUV emission excess (e.g., Bonamente 
et al. 2002, but see also Bowyer et al. 2004 for an alter- 
native explanation); the secondary electrons produced by 
the DM annihilation under the assumption that DM is 
mainly constituted by WIMPs. The overall SZ effect from 
the combination of the different electronic populations in 
Coma can be evaluated properly in the approach described 
by Colafrancesco et al. (2003). The total SZ effect for val- 
ues of M x = 10, 20 and 30 GeV is shown in Fig[U To be 
consistent with the data, the SZ dm effect cannot be due to 
ncutralinos with M x ^ 20 GeV for the assumed value of 
{<tV)a- Since SZdm depends on the quantity {aV)An x , 
the available SZ data set actually an upper limit in the 
(<jV)a — M x plane (see FigEJ • Models with large values 
of (<tV)a and low values of M x which are found in the 
shaded area are excluded by the excess SZdm effect in 
Coma. The limits on {<tV)a set by WMAP restrict fur- 
ther on the available region of the (<jV) a~M x plane. 



3. Conclusions and future prospects 

We have shown that the SZ effect induced by secondary 
electrons produced in XX annihilation is an unavoidable 



consequence of the presence and of the nature of Dark 
Matter in large-scale structures. The analysis of the DM 
induced SZ effect in galaxy clusters provides a comple- 
mentary probe for the presence and for the nature of DM 
in cosmic structures. The available SZ observations on the 
Coma cluster (see Fig|2J) can already set a lower limit to 
the neutralino mass of M x > 17 - 20 GeV (M x > 13 GeV 
at 90 % c.l. with the adopted value of {<tV)a), which 
are consistent with the limits set by accelerators (e.g., 
Belanger et al. 2003). The SZdm signal does not strongly 
depend on the assumed DM density profile at interme- 
diate angular distances from the cluster center and on 
the DM dumpiness since udm is the integral of the total 
Pdm along the line of sight. The presence of a substantial 
SZdm effect is likely to dominate the overall SZ signal at 
frequencies x Si 3.8 — 4.5 providing a negative total SZ 
effect (see Fig^). It is, however, necessary to stress that 
in such frequency range there are other possible contri- 
butions to the SZ effect, like the kinematic effect and the 
non-thermal effect which could provide additional biases 
(see, e.g., Colafrancesco et al. 2003). Nonetheless, the pe- 
culiar spectral shape of the SZdm effect is quite different 
from that of the kinematic SZ effect and of the thermal SZ 
effect and this result allows to disentangle it from the over- 
all SZ signal. An appropriate multifrequency analysis of 
the overall SZ effect based on observations performed on a 
wide spectral range (from the radio to the sub- mm region) 
is required, in principle, to separate the various SZ con- 
tributions and to provide an estimate of the DM induced 
SZ effect. In fact, simultaneous SZ observations at ~ 150 
GHz (where the SZdm deepens the minimum with respect 
to the dominant thermal SZ effect), at ~ 220 GHz (where 
the SZdm dominates the overall SZ effect and produces a 
negative signal instead of the expected ~ null signal) and 
at Si 250 GHz (where the still negative SZ^m decreases 
the overall SZ effect with respect to the dominant thermal 
SZ effect) coupled with X-ray observations which deter- 
mine the gas distribution within the cluster (and hence 
the associated dominant thermal SZ effect) can separate 
the SZ^m from the overall SZ signal, and consequently, 
set constraints on the neutralino mass. Observations of 
the radio-halo emission in the cluster can provide an es- 
timate of the cosmic-ray electron population and conse- 
quently an estimate of the associated non-thermal SZ ef- 
fect (which is usually quite small and with a different spec- 
tral shape at high frequencies, see e.g., Colafrancesco et al. 
2003). The high sensitivity planned for the future SZ ex- 
periments, especially at frequencies x 2.5 and x Si 3.8, 
where the SZdm more clearly manifests itself, can pro- 
vide much stringent limits to the additional SZ effect in- 
duced by DM annihilation. In this context, the next com- 
ing PLANCK-HFI experiment has enough sensitivity to 
probe in details the contributions of various SZ effects in 
the frequency range x ~ 2 — 5. Because the amplitude of 
the SZdm effect increases with decreasing values of M x , 
the high-sensitivity SZ experiments have - hence - the pos- 
sibility to set reliable constraints to the nature, amount 
and spatial distribution of DM in galaxy clusters. We will 
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present elsewhere (Colafrancesco 2004, in preparation) a 
more extended analysis of DM models in the context of 
SZ observations. 

An exciting possibility in this context could be of- 
fered by nearby (with small or zero peculiar velocity) sys- 
tems which are gravitationally dominated by Dark Matter, 
which contain little or no gas (in either hot or warm 
forms) and show absence of non-thermal phenomena con- 
nected with the presence of cosmic rays. In such ideal 
DM systems, the major source of SZ effect would be just 
the one due to the annihilation of the WIMPs. Systems 
which could be assimilable to the ideal "pure" DM halos 
are dwarf spheroidal galaxies and/or low surface bright- 
ness galaxies. These systems seem to be ideal sites for 
studying the DM annihilation indirect signals which re- 
veal themselves in a variety of astrophysical phenomena, 
whose main imprint is the specific gamma-ray emission 
(see, e.g., Evans et al. 2003). In such a context, the possi- 
ble detection of the DM induced SZ effect will provide an 
important complementary approach which can be studied 
by more traditional astronomical techniques. 
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